Quantum Hall effect and the topological number in graphene 
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Recently unusual integer quantum Hall effect was observed in graphene in which the Hall con- 
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ductivity is quantized as Oxy — (±2, ±6, ±10, • • • ) x where e is the electron charge and h is 
the Planck constant. To explain this we consider the energy structure as a function of magnetic 
field (the Hofstadter butterfly diagram) on the honeycomb lattice and the Streda formula for Hall 
conductivity. The quantized Hall conductivity is obtained to be odd integer, ±1, ±3, ±5, • • • times 
two (spin degrees of freedom) when a uniform magnetic field is as high as SOT for example. When 
the system is anisotropic and described by the generalized honeycomb lattice. Hall conductivity can 
be quantized to be any integer number. We also compare the results with those for the square lattice 
under extremely strong magnetic field. 

PACS numbers: 73.43.Lp, 71.70.Di, 81.05.Uw. 



Quantum Hall effect is one of the most spectacular phe- 
nomena in condensed matter physics. The integer quan- 
tum Hall effect has been observed in two-dimensional 
electrons in semiconductor heterostructureQ, and in 
quasi-one-dimensional organic conductors in magnetic- 
field-induced spin density wave stately, Q IE Q- Re- 
cently, quantized Hall effect is observed in graphene 13,1^. 
Unlike the previously studied systems, graphene has a 
unique property that there exist zero modes (so-called 
massless Dirac fermions) in the absence of a magnetic 
field at half-filling due to the honeycomb lattice struc- 
ture. Hall conductivity in graphene has been observed to 
be 




FIG. 1: Generalized honeycomb lattice. Unit vectors are 



vi 



-\a) and V2 



(^a, ^a), where a is the lattice 



constant (a = 0.25nm). A yellow parallelogram is a supercell 
in the presence of a uniform magnetic field of = |(/)o per 
unit cell. 
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where n is an integer. Gusynin and Sharapov|ld| have ex- 
plained these unusual quantized Hall effect as a result of 
four times the quantum Hall conductivity (n + l/2)e^//i 
for each Dirac fermions, where the factor four comes from 
the spin degrees of freedom times number of the Dirac 
fermions in the Brillouin zone. Although their expla- 
nation gives the correct quantum number, their argu- 
ment may be justified only in the low magnetic field limit 
and the logic of quantum Hall conductivity usin^Dirac 
fermions is not correct in general[n|. Peres et al.jE^ ob- 
tained the same result by solving tight-binding electrons 
on a finite size of the honeycomb lattice in a magnetic 
field numerically. 

In this letter we give an approach to the unusual quan- 
tum Hall effect in graphene in terms of topological sta- 
bility. Zero modes at = is a consequence of the 
honeycomb lattice structure which was at best implicit 
in the Dirac Fermion approach. When periodic structure 



is taken into account, there exist gaps in the presence of a 
magnetic field. When chemical potential is in one of the 
gaps, the Hall conductivity is quantized to be an integer, 
which is a topological number and can be obtained by 
Diophantine eg nation |l5 Hsj . The Hall conductivity 
is also obtained by the Streda formula [l^ 
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where N is the total density of states below the gap. Note 
that Hall conductances obtained by the Streda formula 
and by topological numbers are the same since both start 
with the linear response theory. Applying the Streda for- 
mula in the generalized honeycomb lattices in a magnetic 
field, we obtain the unusual Hall effect. This unusual Hall 
effect is also seen in the square lattice at half- filling with a 
magnetic field close to half flux quantum in plaquette^3|- 
The unit cell of the generalized honeycomb lattice con- 
tains two sublattices as shown by the white circles and 
black circles in Fig. ^ There are three non-equivalent 
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FIG. 2: 3D plot of the energy of the honeycomb lattice. Zero 
modes exist at k = ^(7^'^!?) ^ ±(0, |j) 



hopping matrix elements, t^, ^5 and In the usual hon- 
eycomb lattice ta = h = tc^ but they may be different for 
each other under the uniaxial pressure. In the absence of 
the magnetic field zero modes exists, if the condition 



(3) 



IS fulfilledHl. The energy for ta = th = tc is plotted in 
FigEl 

We study a system in a uniform magnetic field with 
a magnetic flux ^ = per unit hexagon, where 

00 = ^nh/e = 4.14 X 10^^ Tm^ is the flux quantum. 
The magnetic field is B = (p/S, where 5* = is the 

area of unit cell of the honeycomb lattice. By choosing a 
suitable gauge and a local gauge transformation, the ef- 
fect of the magnetic field can be treated as a phase factor 
in one of the hopping matrix elements, tc for example, as 
shown in Fig. ^ We take q times larger super cell and 
magnetic Brillouin zone is ^ times smaller in the presence 
of a magnetic field as shown in Figs. ^ and 01 The 
energy is obtained numerically by getting the eigenvalues 
of 2q X 2q matrix for each wave vector k in the magnetic 
Brillouin zone. We take the sufficiently large values of q 
{q = 331 in Figs. 110 and El q = 1201 in Figs. 0111 and El 
and q = 701 in Fig. lTOj) . and then the width of each band 
is very narrow. In this case we do not have to vary k to 
get the energy spectrum as a function of the magnetic 
field. Then taking p = 1, 2, • • • , we get so-called Hof- 
stadter butterfly diagrani,1 9J for generalized honeycomb 
lattice and square lattice[2fl,|2l|. In these figures we also 
plot the marks which shows the p-th, 2p-th, 3p-th, . . ., 
np-th energy bands from the bottom in each = |- We 
also mark the bands of {np + mq)-th. from the bottom, 
where m is positive or negative integer. If there is a gap 
above the {np + mq)-th. band and the chemical potential 




FIG. 3: Brillouin zone for the honeycomb lattice. Gi = 
and G2 = White circles are F points. 

A red hexagon is a first Brillouin zone. A Brillouin zone can 
be taken as the light green diamond. The zero modes occur 
at the corner of the first Brillouin zone (Blue circles) in the 
absence of a magnetic field. A dark green parallelogram is a 
magnetic Brillouin zone in the presence of a uniform magnetic 
field of = 1 00 per unit cell. 
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FIG. 4: Energy spectrum with a magnetic fiux 
eycomb lattice. 



in the hon- 



is in that gap, the Hall conductivity is obtained as 

^2 
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from the Streda formula Eq. Q. It is seen that there 
are gaps corresponding to the Landau levels near the bot- 
tom of the energy at ^ ~ in the generalized honeycomb 
lattices (FigsEUHlandlHI) and square lattice (Fig. Al- 
though each Landau level near the bottom of the energy 
is broadened and splitting into p subband due to the peri- 
odicity of the system in general, the broadening and the 
mini gaps can be neglected in the weak magnetic field 
limit. Therefore, if the chemical potential is fixed near 
the bottom of the band in Fig. usual quantum 
Hall eflFect axy — 2ne^/h with n = ...,3,2,1,0 is ob- 
served as the magnetic field is increased, which can be 
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FIG. 5: The same as Fig.glfor (1)^0 and E ^0. 



Honeycomb t^=1 , \^=^ , 1^=1 .5 




FIG. 6: Energy spectrum with a magnetic flux (j) in the gen- 
erahzed honeycomb lattice with = ^6 = 1, = 1.5 



seen in the lower left part in Figs. ^ and [51 Here we 
have assumed the spin degeneracy and multiply a factor 
of 2. 

Interesting phenomena are seen near half-filling {E ^ 
0). For the usual honeycomb lattice, visible gaps open 
only for the odd number n (Fig.jS}. The band for n = 2 
and 3 are overlapped in the scale of the figures. Although 
there are gaps between every 2q bands when = |(/)o 
in general, only gaps above {np + q)-th. band with odd 
integer n are seen in the E ^ region in that scale. 
Therefore, if the chemical potential is fixed near E = 
(for example jhq in Figs. ^ and [51), quantum Hall effect 
with (Jxy = 2ne^ /h with n = • • • , 7, 5, 3, 1, as the mag- 
netic field is increased. If we change parameter tc in 
the generalized honeycomb lattice, gaps above the even 
quantum number are also visible at finite 0. The gap 
between n = 2 and n = 3 begins to open at (j) ^ 0.018 
for ta = h = 1, tc = 1.5 (Fig.Q). The gap between n = 1 
and n = becomes open for (j) > 0.03 for ta = h = 1, 
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FIG. 7: The same as Fig. jSlin (1)^0 and E ^0. 



Honeycomb 1^=1 , t[j=1 , X^=2 
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FIG. 8: Energy spectrum with a magnetic flux in the hon- 
eycomb lattice of ta = ^6 = 1, = 2 

tc = 1.5. If the chemical potential is in that gap. Hall 
conductivity is zero. In the graphene case, B = 400T 
corresponds to (j) ^ 0.005(/)o 

When ta = tb = 1 and tc = 2, two zero modes in the 
Brillouin zone merge into a confluent mode at = 0_18i]. 
In that case gaps are visible at every integer number n 
dit E ^0 diiid (j)^0 (Fig.O. 

Similar situation as in the honeycomb lattice with 
ta = h = tc in a small magnetic field is seen in the square 
lattice with half flux per unit cell ((/>/(/>o = 1/2). Fig- 
ure ^1 is the Hofstadter butterfly diagram for the square 
lattice. Note the similarity between E ^ 0, (/) ^ 0.5 
in Fig. [ini and £; ^ 0, 0/^o ^ in Fig H In both 
cases visible gaps open only for odd quantum numbers 
(n = 1,3,5,...). For the square lattice with just half 
flux quantum per unit cell, the supercell is twice as large 
as the square unit cell and there exist two zero modes in 
the magnetic Brillouin zone. 

The quantum Hall effect in the square lattice with 
nearly half flux per unit cell and in half- filling case can be 
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where |^| <C 1 and 
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FIG. 9: The same as Fig.|Hlin (/) ^ and E ^ 0. 
square lattice 
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FIG. 10: Energy spectrum with a magnetic flux in the 
square lattice. Hofstadter butterfly 



discussed as follows. When the flux = is applied in 
the square lattice, there are q bands and the contribution 
of the r-th band to the Hall conductivity is given as 
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where tr and tr-i are integers solutions of Diophantine 
eg nation jl5(J 



r = SrQ^ Up, 
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and to = 0. When r bands from the bottom are filled, 
we get 



^ - Va^*) -t — 



= 1 

Consider the case that electrons are near half-filling, 
r 1 , 



P Po 



where \r]\ <C 1. From Eq. we get 
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Then we can conclude from Eq. (|Tn|) that qo should be 
even number (go = 2m where m is integer) . For qo = 2m 
there are zero modes at £^ = 15]. When Eq. ([TT|) is 
satisfied, the r bands are completely filled and the Hall 
conductivity for nearly-half-filled electrons are quantized 
as m times odd integers, because (1 — 2sr)/po cannot be 
an even number. This result is the same as the quantum 
Hall effect on honeycomb lattice in small magnetic field. 

In conclusion we have shown that the unusual quan- 
tum Hall effect in graphene can be explained by the Hof- 
stadter butterfly diagram for the honeycomb lattice. Be- 
cause the honeycomb lattice has two sublattices in a unit 
cell, the similar situation for the half flux quantum per 
unit cell in the square lattice is realized in the small flux 
region. 
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